Microscopic nuclear equation of state with three-body forces and 

neutron star structure 



M. Baldo, G.F. Burgio 
Dipartimento di Fisica, Universitd di Catania and I. N.F.N. Sezione di Catania, c.so Italia 57, 

1-95129 Catania, Italy 

and I. Bombaci 

Dipartimento di Fisica, Universitd di Pisa and I. N.F.N. Sezione di Pisa, Piazza Torricelli 2, 

1-56100 Pisa, Italy 



1 



Abstract 



We calculate static properties of non-rotating neutron stars (NS's) using a 
microscopic equation of state (EOS) for asymmetric nuclear matter. The EOS 
is computed in the framework of the Brueckner-Bethe-Goldstone many-body 
theory. We introduce three-body forces in order to reproduce the correct sat- 
uration point of nuclear matter. A microscopic well behaved EOS is derived. 
We obtain a maximum mass configuration with M^ax = l-SM©, a radius 
R = 9.7 km and a central density nc = 1.34 fm~^. We find the proton frac- 
tion exceeds the critical value x^^'^"', for the onset of direct Urea processes, at 
densities n > 0.45 fm~^. Therefore, in our model, NS's with masses above 
j^Urca ^ O.96M0 can undergo very rapid cooling depending on whether or 
not nucleon superfluidity in the interior of the NS takes place. A comparison 
with other microscopic models for the EOS is done, and neutron star structure 
is calculated for these models too. 



PACS numbers: 97.60.Jd, 21.65.+f 
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In the next few years it is expected that a large amount of novel informations on neu- 
tron stars (NS's) will be available from the new generation of X-ray and 7-ray satellites. 
Therefore, a great interest is devoted presently to the study of NS's and to the prediction of 
their structure on the basis of the properties of dense matter. The equation of state (EOS) 
of NS matter covers a wide density range, from ~ 10 g/ cm^ in the surface to several times 
nuclear matter saturation density (po ~ 2.8 10^^ g/cm^) in the center of the star The 
interior part (core) of a NS is made by asymmetric nuclear matter with a certain lepton 
fraction. At ultra-high density, matter might suffer a transition to other exotic hadronic 
components (like hyperons, a K~ condensate or a deconfined phase of quark matter). The 
possible appearance of such an exotic core has enormous consequences for the neutron star 
and black hole formation mechanism p[. Unfortunately large uncertainities are still present 
in the theoretical treatment of this ultra-dense regime 0,^ . Therefore, in the present work, 
we consider a more conventional picture assuming the NS core is composed only by an 
uncharged mixture of neutrons, protons, electrons and muons in equilibrium with respect 
to the weak interaction (/5-stable matter). Even in this picture, the determination of the 
EOS of asymmetric nuclear matter to describe the core of the NS, remains a formidable 
theoretical problem 0. 

Any "realistic" EOS must satisfy several requirements : i) It must display the correct 
saturation point for symmetric nuclear matter (SNM); ii) it must give a symmetry energy 
compatible with nuclear phenomenology and well behaved at high densities; iii) for SNM 
the incompressibility at saturation must be compatible with phenomenology on monopole 
nuclear oscillations 0]; iv) both for neutron matter (NEM) and SNM the speed of sound 
must not exceed the speed of light (causality condition), at least up to the relevant densities; 
the latter condition is automatically satisfied only in fully relativistic theory. 

In this letter we present results for some NS properties obtained on the basis of a mi- 
croscopic EOS, recently developed which satisfies requirements i-iv, and compare them 
with the predictions of other microscopic EOS's. The Brueckner-Hartree-Fock (BHF) ap- 
proximation for the EOS in SNM, within the continuous choice P], reproduces closely the 
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Brueckner-Bethe-Goldstone (BBG) results which include up to four hole line diagram con- 
tributions ip, as well as the variational calculations |T^, at least up to few times the satu- 
ration density. Non-relativistic calculations, based on purely two-body interactions, fail to 
reproduce the correct saturation point for SNM. This well known deficiency is commonly 
corrected introducing three-body forces (TBF). Unfortunately, it seems not possible to re- 
produce the experimental binding energies of light nuclei and the correct saturation point 
accurately with one simple set of TBF ||10|. Relevant progress has been made in the theory 



of nucleon TBF, but a complete theory is not yet available. In ref. [|TOl a set of simple TBF 
has been introduced within the variational approach. We introduced similar TBF within 
the BHF approach, and we have adjusted the parameters in order to reproduce closely the 
correct saturation point of SNM, since for NS studies this is an essential requirement, and 
there is no reason to believe that TBF be the same as in light nuclei. The corresponding 
EOS (termed BHF3) is depicted in Fig. 1, in comparison with the EOS obtained in BHF 
approximation without three-body forces (BHF2), but using the same two-body force, i.e. 
the Argonne Vn {Avi/^) potential [0). In the same figure, we show the variational EOS 
(WFF) of ref. []IU[ for the 14 + TBF Hamiltonian, and the EOS from a recent Dirac- 



Brueckner calculation (DBHF) |]T^ with the Bonn-A two-body force. The BHF3 EOS 
saturates at Uq = 0.18 fm~'^,E = —15.88 MeV, and is characterized by an incompressibil- 
ity Koo = 240 MeV, very close to the recent phenomenological estimate of ref. 0. In the 
low density region (n < 0.4 /m~^), BHF3 and DBHF equations of state are very similar. At 
higher density, however, the DBHF is stiffer than the BHF3. The discrepancy between these 
two models for the EOS can be easily understood by noticing that the DBHF treatment is 
equivalent |13| to introduce in the non-relativistic BHF2 the three-body force corresponding 
to the excitation of a nucleon-antinucleon pair, the so-called Z-diagram [|T^. The latter is 
repulsive at all densities. In BHF3 treatment, on the contrary, both attractive and repulsive 
three-body forces are introduced and therefore a softer EOS can be expected. 
Fractional polynomial fits to each one of these EOS's allow to compute the corresponding 
pressure and speed of sound Cg to compare with the speed of light c. The ratio Cg/c for all four 



EOS's as a function of the number density is reported in Fig. 2. WFF model (circles) violates 
the causality condition at densities encountered in the core of stars near the maximum mass 
configuration for that model (see fig. 4b) and tab. I). The DBHF calculations need an 



extrapolation to slightly higher densities than the largest one considered in ref. [|T^. The 
extrapolation was done in such a way to keep the causality condition fulfilled. The same 
procedure was followed for BHF3. In the latter case the BHF procedure was well converging 
up to densities n = 0.76 /m'^ for SNM and n = 0.912 /m'^ for NEM. For DBHF the 
causality condition was fulfilled in the extrapolated region only if particular choices of the 
fitting parameters were used, while for the BHF3 the results were insensitive to a wide range 
of variation of the parameters . 



It has to be stressed that the /5-stable matter EOS is strongly dependent on the nuclear 
symmetry energy, which in turns affects the proton concentration |[T^. The latter quantity 



is crucial for the onset of direct Urea processes [0, whose occurrence enhances neutron star 
cooling rates. In our approach, from the difference of the energy per particle E/A in NEM 
and SNM the symmetry energy Egym can be extracted assuming a parabolic dependence on 
the asymmetry parameter /3 = (n„ — Up) / (n„ + Up), being n„ and Up respectively the neutron 
and proton number density. This procedure turns out to be quite reliable [|TB[. The values of 



Egym for the different EOS's are reported in Fig. 3, together with the corresponding proton 
fraction x = {1 — j3)/2. We notice that in both relativistic and non-relativistic Brueckner- 
type calculations, the proton fraction can exceed the "critical" value x^^'^"' = (11 — 15)% 
needed for the occurrence of direct Urea processes ||1^. This is at variance with the WFF 



variational calculation (Fig. 3, circles), which predicts a low absolute value both for the 
simmetry energy and the proton fraction with a slight bend over. For BHF3 model we find 
^Urca _ X3.6%, which correspond to a critical density n^'''^" = 0.447 fm~^. Therefore, BIIF3 
neutron stars with a central density higher than n^'^'^°- develop inner cores in which direct 
Urea processes are allowed. 

The EOS for /^-stable matter can be used in the Tolman-Oppenheimer-Volkoff |T^ 
equations to compute the neutron star mass and radius as a function of the central density. 



For the outer part of the neutron star we have used the equations of state by Feynman- 



Metropohs- Teller |T9[ and Baym-Pethick-Sutherland ^0|], and for the middle- density regime 
(0.001 fm~^ < n < 0.08 fm^^) we use the results of Negele and Vautherin |^T|. In the 
high-density part {n > 0.08 fm~^) we use alternatively the three EOS's discussed above. 
The results are reported in Fig. 4. We display the gravitational mass M^, in units of 
solar mass Mq {Mq = 1.99 10'^'^ g), as a function of the radius R (panel (a)) and the central 
number density Uc (panel (b)). As expected, the stiffest EOS (DBHF) we used in the present 
calculation gives higher maximum mass and lower central density with respect to the non- 
relativistic Brueckner models. The maximum NS mass for the BHF3 is intermediate between 
BHF2 and DBHF, but closer to the latter. The difference between BHF3 and WFF neutron 
stars reflects the discrepancy already noticed for the EOS and mainly for the symmetry 
energy. This point will be discussed in more details in a forthcoming paper ||15|. Table 



I summarizes maximum mass configuration properties, for the different EOS's used in the 
present work. 

In conclusion, we computed some properties of NS's on the basis of a microscopic EOS 
obtained in the framework of BBG many-body theory with two- and three-body nuclear 
interactions. BHF3 EOS satisfies the general physical requirements (points i-iv) discussed 
in the introduction. This is the main feature which distinguishes our BIIF3 EOS with 



respect to other microscopic non-relativistic EOS |10,22|. The calculated maximum mass 



is in agreement with observed NS masses [^]. We found that the neutron star core is 



"proton rich". In fact, the proton fraction at the center of the maximum mass configuration 
in the BHF3 model is x = 26%. Our BHF3 neutron stars with mass above the critical 
value M^''™ = Mdn'^^'^"') = O.96M0 develop inner cores in which direct Urea processes 
can take place. These stars cool very rapidly or not depending on the properties of nuclear 
superluidity (values of the superfluid gaps, critical temperatures, density ranges for the 
superfluid transition) [^,^. Our EOS offers the possibility for a selfconsistent microscopic 
calculation for both the neutron star structure, and nuclear superfluid properties within the 
same many-body approach and with the same nuclear interaction. 
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TABLES 



EOS 




R (km) 




DBHF 


2.063 


10.39 


1.13 


BHF3 


L794 


9.74 


1.34 


BHF2 


1.59 


7.96 


1.95 


WFF 


2.13 


9.4 


1.25 


TABLE L Parameters of the maximum mass configuration: the ratio Mq/Mq 
several EOS's vs. the corresponding radius R and central number density ric- 


is shown for 
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FIGURES 

FIG. 1. The energy per baryon E/A is plotted vs. the number density n for symmetric matter (panel 
(a)) and for neutron matter (panel (b)). Several EOS's are shown, i.e. non-relativistic Brueckner cal- 
culations without (short dashes, BHF2) and with three-body forces (solid line, BHF3) and a relativistic 
Dirac-Brueckner one (long dashes, DBHF). For comparison a variational calculation is also reported (cir- 
cles, WFF). 

FIG. 2. The ratio Cs/c is plotted as a function of the number density for pure neutron matter. The 
different curves refer to BHF2 (short dashes), BHF3 (solid line), DBHF (long dashes) and the variational 
calculations of ref. (circles). 

FIG. 3. The simmetry energy and the proton fraction are shown vs. number density respectively in 
panel (a) and (b). The different curves refer to BHF2 (short dashes), BHF3 (solid line), DBHF (long dashes) 
and the variational calculation WFF (circles). 

FIG. 4. The gravitational mass Mq, expressed in units of solar mass Mq, is displayed vs. radius R and 
the central number density Uc. The notation is the same as in previous figures. 
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E/A (MeV) 




